Scattering of hole excitations in a one-dimensional spinless quantum liquid 
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Luttinger liquid theory accounts for the low energy boson excitations of one-dimensional quantum 
liquids, but disregards the high energy excitations. The most important high energy excitations are 
holes which have infinite lifetime at zero temperature. At finite temperatures they can be scattered 
by thermally excited bosons. We describe the interaction of the hole with the bosons by treating 
it as a mobile impurity in a Luttinger liquid. This approach enables us to evaluate the scattering 
probability at arbitrary interaction strength. In general, the result is expressed in terms of the 
hole spectrum, its dependence on the density and momentum of the fluid, and the parameters of 
the Luttinger liquid Hamiltonian. In the special case of Galilean invariant systems the scattering 
probability is expressed in terms of only the hole spectrum and its dependence on the fluid density. 
We apply our results to the problem of equilibration of one-dimensional quantum liquids. 

PACS numbers: 71.10.Pm 



I. INTRODUCTION 

The Luttinger liquid theory describes low energy 
properties of one-dimensional systems of interacting 
fermionsi^ or bosons^ in terms of excitations with Bosc 
statistics. The latter have an acoustic spectrum and cor- 
respond to long wavelength density fluctuations. De- 
spite the simplicity of the Hamiltonian, the Luttinger 
liquid theory successfully describes the nontrivial power 
law behavior 5 of various correlation functions of one- 
dimensional systems at small frequencies and wavenum- 
bers. Many of them, such as the power law energy de- 
pendence of the tunneling density of states, have been 
experimentally confirmed £ 

The power law behavior of correlators in one- 
dimensional systems is not limited to small frequencies 
and wavenumbers. For example, even at large wave num- 
bers q ^ no, where no is the fluid density, the zero tem- 
perature dynamic structure factor (the Fourier transform 
of the density-density correlator) exhibits a power law 
singularity S(q,uj) oc [fun — e(q)] a 9(fno — e(q)). Its loca- 
tion e(q) defines the spectral edge below which the system 
cannot absorb excitations. It corresponds to the lowest 
energy state of the system with momentum Kq. The po- 
sition of the spectral edge e(q) is periodic in q with the 
period 2nno (in the spinless case considered here), as il- 
lustrated in Fig. [TJ 

The nature of states corresponding to the spectral edge 
and periodicity of the latter may be illustrated by the 
simple example of noninteracting fermions in one dimen- 
sion. Because the energy cost associated with trans- 
ferring a particle between the two Fermi points van- 
ishes in the thermodynamic limit, the minimal energy 
states with momenta differing by an integer multiple of 
2hkF = 2nhno are degenerate, resulting in the period- 
icity of the spectral edge. It is thus sufficient to con- 
sider states with wave vectors in the fundamental domain 
< q < 27m . The states at the spectral edge corre- 
spond to hole excitations which are obtained by moving 



FIG. 1: At a given momentum hq the energy of a one di- 
mensional system is limited from below by the spectral edge 
e(q). Thus the imaginary part of the structure factor S(q,u), 
describing the dissipation, vanishes outside the shaded region. 



a fcrmion with the wave vector fcp — q to the right Fermi 
point k — hp. 

In the presence of interactions between the fermions 
the energy cost associated with the transfer of a particle 
between the opposite Fermi points still vanishes in the 
thermodynamic limit, and the spectral edge e(q) remains 
periodic in q with the period Itxtiq. A state at the spectral 
edge with a wave vector in the range < q < 2-ktiq 
may again be viewed as the q = ground state with 
an additional hole. However, the hole is now dressed 
by the interactions, and its energy e{q) is renormalized. 
The above picture of the spectral edge also applies to 
bosonic fluids, where holes are known as Lieb's type-II 
excitations 2- 

At q — > and q — > 2irno the energy e(q) is small and 
the hole can be decomposed into the bosonic excitations 
of the Luttinger liquid. In contrast, the high energy holes 
with q ~ no are not described by the Luttinger liquid the- 
ory and should be viewed as distinct from the Luttinger 
liquid bosons. As noted above, the states with a single 
hole are lowest energy states of the system with a given 
momentum. Therefore they cannot decay into bosons 
and their lifetime is infinite. 

High energy hole excitations may be produced in the 
system by external probes (e.g., optically or by tunnel- 
ing) or by thermal fluctuations. Therefore the problem of 
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FIG. 2: Scattering of a hole at wavenumber Q to a state 
Q + <5Q involves absorption of one bosonic excitations and 
emission of another. 



their dynamics is of considerable interest. Although the 
lifetime of the hole states is infinite at zero temperature, 
the situation changes at T > 0. In this case the hole can 
scatter off thermally excited bosons. Since the resulting 
energy change Se < T, a hole with energy e(q) 3> T re- 
mains distinct from the Luttinger liquid bosons and may 
be treated as a mobile impurity in a Luttinger liquid 4r— 

In this paper we evaluate the probability of scatter- 
ing of a hole by the bosonic excitations. The scattering 
event will be assumed to change the wavenumber of the 
hole Q by a small amount SQ. In order to conserve both 
the energy and momentum of the system, the hole has to 
absorb one boson and emit another, see Fig. [2] Such pro- 
cesses were first considered by Castro Neto and Fisher— 
in the context of the dynamics of mobile impurities in the 
Luttinger liquid. The scattering probability near the top 
of the spectrum, Q — 7m , determines the rate of equi- 
libration of one-dimensional quantum liquids studied in 
Ref. U3- In the limit of strong repulsion the scattering 
processes of Fig. [2] were studied in the context of equili- 
bration of the one-dimensional Wigner crystal , 18 ' 19 while 
for weakly interacting bosons they are responsible for the 
decay of the so-called dark solitonJ^ The mobile impurity 
approach enables us to treat this problem at any inter- 
action strength. We express the scattering probability in 
terms of the hole spectrum e(q) and its dependence on 
the density and velocity of the fluid. For Galilean invari- 
ant systems the result can be expressed in terms of only 
the spectrum and its dependence of the fluid density. 

The standard bosonization of Luttinger liquids ji* 2 - used 
in the previous treatments of mobile impurities ^ 14 ' 20 " — 
corresponds to the Eulerian description of the liquid. In 
this approach the dynamical degrees of freedom describ- 
ing an element of the fluid are given as functions of time 
t and the instantaneous position x of the element in 
the laboratory frame. Alternatively, the liquid may be 
described using Lagrangian variables. In the latter ap- 
proach the dynamical degrees of freedom are labeled by 
the position y of the fluid element in a reference state of 
uniform density uq. For a liquid moving with a uniform 
velocity the transformation from Eulerian to Lagrangian 
variables is equivalent to a Galilean transformation from 
a laboratory frame to a reference frame moving with the 



fluid. Therefore one expects that Galilean invariant sys- 
tems are more naturally described in Lagrangian vari- 
ables. This expectation is borne out: the use of La- 
grangian variables considerably simplifies evaluation of 
the scattering amplitude. 

On the other hand, Galilean invariance is not a uni- 
versal property of Luttinger liquids. For instance, elec- 
trons in solids move in the periodic potential of the lattice 
which gives rise to the band structure of the spectrum. 
In this case the quadratic spectrum p 2 /2m required for 
Galilean invariance appears only near the band edges. 
A similar situation arises when the concept of Luttinger 
liquid is applied to spin chains and one-dimensional sys- 
tems of cold atoms in optical lattices. In the absence of 
Galilean invariance Lagrangian variables offer no obvi- 
ous advantages. To address this regime we develop the 
theory of hole scattering using the standard Luttinger 
liquid theory^ based on the Eulerian approach to fluid 
dynamics. 

The paper is organized as follows. In Sec. [TT] we con- 
sider the Galilean invariant case, and develop theoretical 
description of quantum liquids in Lagrangian variables. 
We obtain the scattering probability of the hole in terms 
of its spectrum. In Sec. IIIII we consider the general case, 
where Galilean invariance is not assumed, and use the 
conventional approach based on the Eulerian variables 
to obtain the scattering probability of the excitation. In 
Sec. II VI we verify that the two approaches give the same 
result in the case of a Galilean invariant system and adapt 
our calculations to the problem of a massive mobile im- 
purity in a Luttinger liquid. We also discuss the impli- 
cations of our results for the problem of equilibration of 
the Luttinger liquid. 



II. SCATTERING OF HOLES IN GALILEAN 
INVARIANT SYSTEMS 

In this section we consider systems which possess 
Galilean invariance. Evaluation of the hole scattering 
probability is simplified by developing the theory in La- 
grangian variables, 26 which are frequently used to de- 
scribe one-dimensional flows in classical hydrodynamics. 



A. Hamiltonian of a quantum liquid in Lagrangian 
variables 

We first consider a liquid in which only low energy 
excitations are present. These excitations are essentially 
sound waves with small wavenumbers q <C uq. In this 
regime the discreteness of particles is not important and 
the liquid may be described as a continuum in the spirit 
of hydrodynamics. 

The Hamiltonian describing the long wavelength exci- 
tations may be obtained by coarse-graining the fluid into 
small elements. The Lagrangian coordinate y of the fluid 
element corresponds to its position in a reference state 
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of uniform density no- I n the presence of sound waves 
the particle density n(y) inside each element may devi- 
ate from no- The length of each element in the reference 
state Ay is assumed to be sufficiently small, qAy -C 1, 
so that the variations of the fluid density inside it are 
negligible. On the other hand, the number of particles in 
the element is assumed to be large, no Ay 3> 1. 

Taking advantage of the Galilean invariance we write 
the energy of the fluid element as a sum of the kinetic 
energy of its center of mass motion and the internal en- 
ergy, 



AE 



(AP) 2 
2mnoAy 



U(n) n Ay. 



(1) 



Here AP is the momentum of the element and m is the 
mass of the particles. In the absence of high energy ex- 
citations the internal energy is given by the energy per 
particle in the ground state U(n) multiplied by the num- 
ber of particles noAy. 

The Hamiltonian of the liquid is written in terms of 
two dynamical variables, the displacement of the fluid 
elements from their reference positions u(y) and the con- 
jugate momentum density p(y) = AP/Ay, satisfying the 
commutation relation [u(y),p(y / )] — ihS(y — y'). Sum- 
ming the energies ([1} of the fluid elements we obtain 



Hr 



Imno 



n U(n) 



dy. 



(2) 



The Lagrangian variable u(y) enters the Hamiltonian via 
the particle density 



n(y) 



no 



l + u'(yY 



(3) 



where prime denotes the derivative. Equation ([3]) follows 
immediately from the relation 



x = y + u{y) 



(4) 



between the physical coordinate x of the fluid element 
and the Lagrangian variable y. 

For the subsequent discussion it is sufficient to expand 
the Hamiltonian of the liquid to third order in the defor- 
mation it', 



Hr 



(J: 



mnov , 2 



\2mno 



u 12 - au ri dy. (5) 



The speed of sound in the liquid v and the coefficient a 
in the anharmonic term are given by 



1/2 



v(n ) = |- [2noU'(no) + n 2 U"(no)]j , (6) 
o(no) - n 2 ) U'(n )+nlU"(n ) + ±nl ) U"'(n ), (7) 

where derivatives are again denoted by prime. 

The quadratic part of Eq. (|S|) is the Luttinger liquid 
Hamiltonian in Lagrangian variables. It can be brought 



to the diagonal form ^ hv\q\b q b q by introducing the bo- 
son operators b q via the standard procedure^ 



2mnoLv\q\ 



(b q e™ + b\e*») , ( 8 ) 



p(y) 



. E y^M {bgeiqy _ ble - iqy) . (9) 



where L is the system size. 

Because of the continuum approximation made in the 
description of the liquid the sums in the above equa- 
tions include only small wavevectors \q\ <C qo, where 
qa = 1 /Ay. Therefore if a high energy hole is present 
in the liquid, it needs to be treated as an additional en- 
tity. 



B. Description of the hole in Lagrangian variables 

Let us now consider the fluid with a single high energy 
hole excitation. The Lagrangian coordinate Y of the hole 
refers to the fluid element containing it. As discussed in 
the introduction the hole excitation is obtained by mov- 
ing a particle from state Uf ~ q to the Fermi point fcp . 
Thus the presence of the hole has no effect on the mass 
of the fluid element, and the expression for its kinetic en- 
ergy given by the first term in Eq. (fTJ) remains unchanged. 
On the other hand, the internal energy increases by the 
excitation energy e(q, n). (Here we take into account the 
dependence of the latter on the fluid density n.) 

It is important to note that in the presence of sound 
waves the physical size of the fluid element Aa; differs 
from its size in the reference state, Aa; — (no/n)Ay. 
As a result the physical wavenumber q differs from the 
wavenumber Q corresponding to the Lagrangian coor- 
dinate of the excitation, q = (n/no)Q- Therefore it is 
convenient to introduce the excitation spectrum e(Q,n) 
with respect to the Lagrangian wavenumber Q, which is 
related to the physical spectrum e(q, n) by 



e(Q, n) = e(nQ/n Q , n). 



(10) 



We now conclude that the presence of a hole at point Y 
is accounted for by the contribution to the Hamiltonian 
in the form 



H h =e{-idY,n(Y)). 



(11) 



The density n(Y) at the location of the hole is affected by 
the boson excitations, see Eqs. ([3]) and ©. This results 
in the interaction of the hole with the bosons. 



C. Scattering probability 

We now consider scattering processes shown in Fig. [5] 
and evaluate their rate Wq.q+sq- To this end we expand 
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the excitation energy to second order in u' 
e {2) (Q,n(Y)) = eg - n d n e Q u'(Y) 



n d n £ Q + -n d n e Q 



u'\Y). (12) 



In this equation and the subsequent results the hole en- 
ergy cq = e(Q, n) and its partial derivatives are evalu- 

(2) 

ated at n = uq. The corresponding Hamiltonian Hi 
is obtained by substituting Q — —idy and symmetrizing 
the operators. Because of the relation §8§ contains 
both linear and quadratic coupling between the bosons 
and the hole. 

The scattering process illustrated in Fig. [5] involves two 
bosons. Its amplitude may be obtained using perturba- 
tion theory in the coupling between the hole and the 
bosons. One contribution to the scattering amplitude 
arises from the quadratic coupling in the second line of 
Eq. (Tj"2"|) taken in the first order of perturbation theory. 
Keeping in mind Eq. ([5]) one easily concludes that such 
amplitude is proportional to the small momentum trans- 
fer SQ. 

Another contribution arises from the linear coupling 
term which scales as \/SQ. The process shown in Fig. [2] 
is realized in the second order perturbation theory in the 
linear coupling. Since the energy denominator is propor- 
tional to SQ one may naively expect the respective am- 
plitude to scale as (\/<^3) /SQ = (SQ) a . However, there 
are two processes in which the two participating bosons 
are created and destroyed in opposite order. Since the 
corresponding energy denominators have opposite signs, 
the amplitudes of such processes cancel each other in 
leading order in SQ. The remaining subleading contri- 
bution again scales as SQ. 

Finally there is a contribution in the second order per- 
turbation theory, which arises from the combination of 
the linear coupling of the bosons to the hole and the au' 3 
term in Eq. (JSJ) - The latter scales as (SQ) 3 ^ 2 , resulting 
in the scattering amplitude linear in SQ. 

One can avoid the somewhat tedious calculation out- 
lined above by performing a unitary transformation 
U^(Hl + Hfp)U of the Hamiltonian, which eliminates 
the linear coupling between the hole and the bosons to 
leading order in SQ. The operator U should be chosen in 
the form 



U = exp l -f u u{Y) + l -f p 



p(y)dyj 



Here the coefficients /„ and f p are given by 
f _ n VQd n €Q _ d n e Q 

■/« — -.2 ..2 ' JP ~ 



V — V, 



m(v 2 — Vq) ' 



(13) 



(14) 



where vq is shorthand notation for the velocity of the 
hole, 



vq = v(Q, n ) = j i dQ£(Q, n). 



(15) 



Because in the above unitary transformation Q denotes 
the initial momentum of the hole rather than the oper- 
ator — idy the linear coupling between the hole and the 
bosons is removed only to leading order in SQ. However, 
in its absence the subleading contribution, which scales 
as (SQ) 3 / 2 , gives a negligible correction to the scatter- 
ing amplitude, [(SQ) 3/2 ] 2 / SQ = (SQ) 2 . The dominant 
contribution to the scattering amplitude arises from the 
terms in the transformed Hamiltonian that couple the 
hole to second powers of boson fields, u' 2 (Y) and p 2 (Y). 
This coupling has the form 

y - 3a/ p + n d n e Q + \^d 2 n <L Q - v Q f u - n d n v Q f u 
^u' 2 (Y)-SrP 2 (Y), (16) 



2m 



2m 



Q 



where the momentum dependent effective mass m*q is 
defined in terms of the curvature of the spectrum 



1 



1 



m* Q h 2 



In addition we introduced the notation 

V Q = v(Q,n) = v(nQ/n ,n) 



(17) 



(18) 



for the velocity of the hole as a function of the Lagrangian 
wavenumber, cf. Eq. (|10[) . Although vq — vq at n — no, 
their dependences on density are different. In particular, 



d n VQ = d n VQ 



HQ 

m* n n Q 



(19) 



at n = no. 

The matrix element for the scattering process in which 
a boson with momentum q\ is annihilated and a boson 
with momentum qi is created is obtained by expressing 
the deformation v! and the momentum density p in the 
above expression in terms of the boson creation and an- 
nihilation operators using Eqs. © and ©, 



mnoLv 
+n d n v Q f., 



3a/ p - n d n e Q - -riQ^eQ + v Q f u 



fu - (mn vf p ) 2 
2m* n 



(20) 



Here we assumed that q\q2 < 0, see Fig. [2j 

The coefficient a may be expressed in terms of the 
density dependence of sound velocity v with the aid of 
Eqs. © and 0, 



3a(n) = mn [ v 2 + —nd n v 2 



(21) 



Using this relation and Eq. (fT"4|) we can express the scat- 
tering amplitude in terms of the spectrum of the hole, 



_ frn ^\qiq 2 \ 
qiq2 ~ 2mLv(v 2 -v 2 Q ) C 



(22) 
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where wc introduced the notation 



The field 4> is defined in terms of the particle density 



T Q = (d n e Q ) d n (v 2 - Vq) ~ (v 2 - v 2 Q ) d 2 n e Q 



(23) 



The scattering probability per unit time is given by the 
Fermi golden rule, 

Wq,q+8 Q = K^\ 2 N qi (N q2 + -12- SQ) 



x<5(e<2 - £Q+«5Q + hv\qi\ - hv\q 2 \), 



(24) 



where N q is the occupation number of the boson state q. 
We have normalized the probability Wq } q+sq in such a 
way that the total scattering rate is given by the integral 

Wq = J W Q , Q+SQ dSQ. (25) 

Using Eq. (|2"2"|) we immediately obtain 



<rSQ 



n 2 a N qi (N q2 +l) 2 



64irm 2 v 5 (v 2 — Vq) 
Here q\ and q 2 are given by the relations 



T$ (5QY 



(26) 



qi = ^SQ+^\5Q\, q2 = -lsQ + ^-\5Q\, (27) 
2 Zv Z Zv 

arising from conservation of energy and momentum. 

Equations (|23|) and (f26|) express the scattering rate 
of a high energy hole excitation in a Galilean invariant 
quantum fluid in terms of the hole spectrum (|10j) . They 
are the main result of this section. 



III. SCATTERING OF HOLES IN THE 
ABSENCE OF GALILEAN INVARIANCE 

We now turn to the more general situation in which 
the quantum liquid is not assumed to possess Galilean 
invariance. As discussed above, in this case Lagrangian 
variables offer no obvious advantages. Therefore we ap- 
ply the standard (Eulerian) theory^ of the Luttinger 
liquid. 



n(x) = no H S7(j)(x) 



(29) 



The field 9 accounts for the motion of the liquid. The 
latter may be characterized by the momentum per par- 
ticle 



k(x) 



1 ( 1 ^ I s 1 
■p(x) +p(x) 



n(x) 



(30) 



2 \n(x) 
Here the momentum density 

= \ ^2ipi s ( x - x i) + s ( x - x i) pi] ( 31 ) 



is defined in terms of the coordinates xi and momenta 
Pi of the physical particles. Using the commutators 
[xi,pi'] = ihdw and definition of the particle density 



(32) 



one easily obtains the commutation relation 

[n(x), k{x')} = -ihV5{x - x'). (33) 

Comparing this relation with Eqs. (|2"5)) and (|2^|) we iden- 
tify 



k(x) = -HV9{x). 



(34) 



Thus the gradient of the boson field 8 determines mo- 
mentum per particle in the fluid. 

It is worth mentioning that in the Galilean invariant 
case k(x) — mV(x) and V#(x) gives the expression 



V(x) = V6(x) 



(35) 



for the fluid velocity. The relation between the fields 
4>{x), 0(x) and the Lagrangian variables u(y), p(y) used 
in Sec. |H] is given by Eqs. ©, O, ([3D]), and as 
well as Eq. Q which expresses the physical coordinate x 
in terms of the Lagrangian coordinate y. 



A. Eulerian description of one-dimensional 
quantum liquids 

The standard theory of a Luttinger liquid describes the 
system by two bosonic fields <fi(x) and 0(x) satisfying the 
canonical commutation relations 

[<j){x), V0(x')] = i^5{x - x'). (28) 

In contrast to the coordinate y in the Lagrangian ap- 
proach of Sec. HH the coordinate x denotes the real space 
position, see Eq. (Q}. 



B. Hamiltonian of the liquid in the presence of the 

hole 

The standard form of the Hamiltonian of a Luttinger 
liquid is 

Ha = z^J dx [^ W ) 2 + ^W) 2 ] ( 36 ) 

(see, e.g., Sec. 3.1 of Ref. 0). Here the velocity v and the 
dimensionless Luttinger liquid parameter K depend on 
the density of the particles and the interactions between 
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them. This Hamiltonian describes noninteracting bosons 
with an acoustic spectrum. It can be brought to the form 



with the help of the relations 



2L 



(37) 



n(g)(6 9 e^-6t e -^),(38) 



(39) 



which express the fields tf> and 6 in terms of the boson 
creation and annihilation operators b q and b q . 

The Hamiltonian Hq describes the low energy proper- 
ties of the system and may be viewed as the fixed point 
Hamiltonian in the renormalization group sense. The 
leading irrelevant perturbation is given by cubic in the 
boson fields corrections r^l 



dx[a e {V<j)){V6) 2 + a^V^f 



(40) 



Here we assumed that the system is invariant under in- 
version, x — > —x. In this case the Hamiltonian cannot 
contain odd powers of V(9. The parameters ctg and 
can be expressed in terms of density dependent param- 
eters of the fixed point Hamiltonian (|36]l. Indeed, ac- 
cording to Eq. (|21?1) a small change of density 5n shifts 
— > Wtfi + nSn. As a result the cubic perturbation (|40)) 
generates a small correction to the quadratic Hamilto- 
nian (|36[) . The corresponding change in its parameters is 
5(hvK /2ir) = TragSn, and 8(hv /2nK) — 3a^5n. We thus 
find the expressions 



a s = d n (vK) . 
2n z 



67T 25n (x) 



(41) 



for the parameters of the Hamiltonian (|4T))) . 

Let us now consider a hole excitation at point X. In 
general its energy e depends not only on the wave number 
Q = —idx but also on the density of the liquid n(X) and 
its motion quantified by k{X). If the spectrum of the hole 
in a uniform liquid s(Q, n, k) is known, the correction to 
the Hamiltonian of the system due to the presence of the 
hole may be written as 



H h = e{-idx,n{X),K{X)). 



(42) 



The right hand side here is assumed to be symmetrized 
with respect to the operators —idx, n(X), and n(X). 
In view of the relations (|2"5|) and (|3~4"|) this Hamiltonian 
describes the hole interacting with the boson fields 4> and 



C. Scattering probability 



obtained in Lagrangian variables for the Galilean invari- 
ant case, Eqs. ([5]) and (fTTj) . The main difference is that 
the spectrum of the hole in the Eulerian description de- 
pends not only on the fluid density n but also on its mo- 
mentum k. The rate Wq^q+sq of scattering of the hole 
by bosons, Fig. [2J can be found by repeating the steps 
of Sec. Ill CI First we perform a unitary transformation 
U^HU of the Hamiltonian with the operator 



U = exp [if^iX) + ifg0(X)} . 



(43) 



Upon such a transformation the Hamiltonian takes the 
form Ho + H a + Hh where 

H h = e(-id x + iW + fgV6,n + V<b/n,-hV8) 



+hvKf <p V9+—fgV<b 
K 

+27r/ a e V<£V0. 



(44) 



Here the gradients of the boson fields are evaluated at 
point X. 

The parameters and fg should be chosen such that 
the linear terms in the expansion of Eq. (|44[) in powers of 
V(j) and V9 vanish to leading order in — idx — Q, where Q 
is the initial wavenumber of the hole, Fig. [21 This yields 



f'b 



1 VQKd n SQ + irhv d K e Q 

V z — Vn 



fe = - 



TThK 

1 vKd n £ Q + irhvQd K EQ 



V 2 — V 2 



(45) 
(46) 



Q 



where we used the shorthand notations eq — e(Q,n,n) 
and vq = dQEQ/h. In Eqs. (|45|). (|46|). and the subsequent 
results Eq and its partial derivatives are evaluated at n — 
no and k = 0. 

With the above choice of and fg the Hamiltonian 
(1441) takes the form 



h 2 



H h = 




2ir 2 7T 



9 n VQ 



(47) 



Here we omitted linear in V0 and V# terms with coeffi- 
cients small in —idx~Q, the quadratic term proportional 
to (V0)(V0), and higher powers of Vcj) and V8. Such 
perturbations do not affect the amplitude of the scatter- 
ing process depicted in Fig. [2] to leading order in SQ. In 
Eq. (|47|) the effective mass of the hole mJ, is defined by 



1 



1 



<9q£(<3,"o,0). 



(48) 



The Hamiltonian of the liquid in the presence of the 
hole given by Eqs. (|3l)]). (I4U1) and (|4"2"j) is similar to that 



We now substitute relations (J3S]) and (J2HJ) into Eq. (|4"7) 
and extract the matrix element t qiq2 corresponding to the 
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scattering process in Fig. [2j Using Eqs. (l4"T1) . (l45j) . and 
(|46[) we express the matrix element i gig2 in the form 



K V|gig 2 | 

2ttL V 2 - V% 



(49) 



where 



1 



71" /l 



-2 ( ^d K v Q - -^9„i^Vw9„£Q + v Q ^rd K e Q 



-2(d n v Q ) 



VQd n e Q +v^rd K e Q 



(50) 



Substituting Eq. (|49|) into the Fermi golden rule expres- 
sion (|24|) we obtain the scattering rate 



W< 



K 2 N qi (N q2 + l) 



(51) 



where g x and are again given by Eq. (|2"T|) . 

Equations ([ST]) and (|50|) express the scattering rate 
of the hole in terms of its spectrum e(Q, n, k) and the 
parameters of the Luttinger liquid. They are the main 
result of this section. 



IV. DISCUSSION OF THE RESULTS 

In sections QT] and IIIII we obtained expressions for 
the rate Wq^q+sq, of scattering processes illustrated in 
Fig. [2] In Sec. [TT] we considered Galilean invariant sys- 
tems and obtained the hole scattering rate in the form 
of Eqs. (|26|) and (|23|) . Our consideration in Sec. IIIII did 
not assume Galilean invariance and yielded the scatter- 
ing rate in a somewhat more complicated form, Eqs. (|5"Tj) 
and ([50)) . Below we compare these results and discuss 
their applications to the problem of equilibration of one- 
dimensional quantum liquids and to the problem of dy- 
namics of a mobile impurity in a Luttinger liquid. 



Galilean invariant systems 



To obtain the expressions (f2l)]) and (|23p for the scat- 
tering rate we described the liquid in Lagrangian coor- 
dinates. The state of the hole was parameterized by the 
Lagrangian wave number conjugate to the Lagrangian 
coordinate of the hole Y, rather than the conventional 
wave number conjugate to the spatial coordinate X. Ac- 
cordingly, the excitation energy eg — e{Q, n) and velocity 
vq = v(Q,n) were expressed in terms of the Lagrangian 
wavenumber using Eqs. (|10|) and (JTJ 



To illustrate this point let us consider the simple spe- 
cial case of noninteracting fermions. Recalling that the 
hole excitation is created by moving a fermion from state 
kp — Q to state kp we find its energy and velocity in the 
form 

h 2 h 
SQ = ^Q(2nn-Q), w Q = -(tto - Q). (52) 
2m m 

To express the energy and velocity in terms of the La- 
grangian wave number we substitute Q — > hQ/uq, see 
Eqs. (HOJ) and (HI]), and find 



2mriQ 



Q(2irn - Q), v Q 



Tin 
mno 



(tttio-Q). (53) 



Substituting (|53]l into f[23]) we obtain Tq = 0. This is the 
expected result, as no scattering of excitations may occur 
in a system of noninteracting fermions. An erroneous 
substitution of the eq and vq given by Eq. (|52|) for cq 
and vq into Eq. (|2"3"|) would result in Tq ^ 0. 

Although less common than the Eulerian variables in 
hydrodynamics, Lagrangian variables arise naturally in 
the theory of elasticity. In particular, in the case of 
a one-dimensional Wigner crystal (i.e. quantum anhar- 
monic chain) the Lagrangian coordinate is essentially the 
number of a site of the Wigner lattice. Thus the phonon 
spectrum loq of the Wigner crystal is naturally expressed 
in terms of the Lagrangian wavenumber. The weakness of 
interactions between the phonons in the Wigner crystal 
enables one to develop a microscopic theory of scatter- 
ing of a high energy phonon by acoustic modes i 18 ! 19 On 
the other hand, the Wigner crystal is simply the limiting 
case of a quantum fluid in the regime of extremely strong 
repulsion between the particles. Thus the scattering of a 
high energy phonon can also be studied phenomenologi- 
cally using the approach of Sec. |TTJ Indeed, substituting 
the phonon energy e(Q,n) = hujQ into Eq. (|23[) we re- 
cover the result of Ref. [jjj 

Experimentally, it is easier to measure the excita- 
tion spectrum e(Q, n) as a function of the conventional 
wavenumber. To obtain the scattering rate in terms of 
e(Q, n) using Eq. (l26l) one should substitute the relations 
{TDJ and CEH]) into Eq. ([23]). This yields, 



(d n e Q ) d n {v 2 - v Q ) - (v 2 - v Q ) dls Q 



nr 



{dnSQ? 



2hvQ 
n 



(v Q d n v - vd n v Q ) . 



(54) 



As before, the right hand side here is evaluated at n = ng. 
Substitution of the expressions (f52"j) for the energy and 
the velocity of the hole in a noninteracting Fermi gas into 
Eq. (|5i|) gives the correct result Tq = 0. 

Our discussion of hole scattering in Sec. IIIII did not 
assume Galilean invariance of the system. The resulting 
scattering rate Wq^q+sq is given by Eqs. (l5"Tj) and ([501) . 
For systems that possess Galilean invariance these results 
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should agree with those of Sec. [n] Eqs. and (l23l) . 
The Luttinger liquid parameter K in Galilean invariant 
systems can be expressed in terms of the particle density 
n and the velocity v of excitations as 



K = 



it Yin 



(55) 



Substituting this expression into Eq. (|5Tj) and setting 
n = n we recover Eq. ([2"6"]) . To demonstrate that the 
resulting scattering rates are equal, we also need to de- 
rive the expression ([23]) for Tq from Eq. ([50]) . 

In Sec. Mil the hole was described by the dependence 
of its energy on the wavenumber Q, particle density n 
and the momentum per particle n. In Galilean invariant 
systems k = mV, where V is the fluid velocity. The 
energy of the excitation e(Q, n, k) in the moving fluid 
differs from its energy e(Q,n) in the stationary fluid by 
HQV (see e.g. Ref. [28i), resulting in 



e(Q, n, k) — e(Q, n) + HQ- 



(56) 



Substituting Eqs. ([55]) and (151))) into Eq. ([5TJ]) we obtain 
the result (|54|) . which is equivalent to Eq. ([23]) . 



B. Hole vs. particle-hole excitation 

Throughout this paper we considered hole excitations 
created by moving a particle from a state kp — Q to 
the Fermi point kp. On the other hand, our discussion 
of hole scattering in Sec. IHII did not rely on this physical 
picture and thus should be applicable to a hole excitation 
obtained by removing a particle from the system. The 
wavenumber Q of such an excitation is related to Q by 



Q(n, k) = Q— k R (n, n), 



(57) 



where k R is the wavenumber of the particle at the right 
Fermi point. Similarly, its energy i is related to e by 



e(Q, n, k) = /J. R (n, n) + e(Q — k R {n, n),n, n), 



(58) 



where ^l r is the energy of the particle at the right Fermi 
point. Since the two physical pictures of the hole excita- 
tion are equivalent, one should expect to find the same 
scattering rate for the particle-hole excitation with en- 
ergy e(Q, n, k) as for the hole (missing particle) with en- 
ergy e(Q - im ,n, k). 

To verify such a feature of Eq. ([50]) we need to ob- 
tain the expressions for k R (n,K) and the derivatives of 
/j, r (ti,k). To this end we express n and k in terms of 
the numbers N R and N L of the right- and left-moving 
particles in a uniform system, 



N R + N L „ N R - N L 

n= , k = nn . 

L L 

Then from k F = (2tt/L)N r we find 



(59) 



To find the derivatives of /i r (ti,k) we substitute V(f> = 
7r(n — no) and V# = —k/H into the Hamiltonian of the 
liquid given by Eqs. (|36|) and ([40]) . Differentiating the 
resulting expression for the energy of uniform liquid with 
respect to N R we find the expression 

[i R {n,K) — [i R (no, 0) + ^r-{n - n Q ) + vKn 



K 



nag , 2 w 
— 5- Ik + ZTrnin — no)K\ 
nr 

+3n 3 a^(n — n ) 2 + . . . 



(61) 



valid to second order in n — no and k. From the first line 
one then immediately obtains the first derivatives of n R 
in the form 



it 



TTh—, d K IJL R 
A 



vK. 



(62) 



The most convenient expression for the second deriva- 
tive dn/j, R is obtained simply by differentiating the first 
of the expressions (|62p . To find d 2 /i R we notice that ac- 
cording to the second line of Eq. (fBTj) we have d 2 /i R = 
{nh)- x d n {d K iJL R ). Then from Eq. (jMJ) we obtain 



d„u K = Trhd„ 



2,,-R 



A" 



d n (vK) 

7th 



(63) 



Using the expressions (|60|) . ([62]). and ([63]) one can show 
that Tq obtained by the substitution of e(Q, n, k) in the 
form ([5"5]) into Eq. ([50]) coincides with that obtained by 
simply replacing e{Q, n, k) — > e{Q—irnQ, n, k) in Eq. (|50|) . 



C. Scattering of a mobile impurity in a Luttinger 
liquid 

Apart from a hole excitation created by either mov- 
ing a particle to the Fermi point or removing it from the 
system, the results of Sec. IIIII apply to any mobile impu- 
rity in a Luttinger liquid. Our result for the scattering 
rate given by Eqs. ([5T]) and ([5TJ]) is applicable regard- 
less of Galilean invariance of the system, but requires the 
knowledge of the energy of the mobile impurity e(Q, n, k) 
as a function of the density and momentum of the liq- 
uid. Similar to the case of a hole excitation, Eq. (|56p . 
this relation simplifies in the presence of Galilean invari- 
ance. One should note, however, that unlike an intrinsic 
particle-hole excitation, a foreign particle has a non-zero 
mass M, which affects the dynamics of the liquid. The 
dependence of the energy s p (V) of the impurity with mo- 
mentum p on the velocity V of the fluid can be obtained 
by comparing the expressions for the energy of the sys- 
tem in the stationary frame and that moving with the 
fluid, 



s p (V) = e p -Mv(Q)+pV 



-MV 



(64) 



k R (n, k) — nn + 



(60) 



see Ref. 2i 
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Denoting the momentum of the impurity p = HQ, the 
velocity of the fluid V — k/tu, and expanding to second 
order in k, we obtain 



s(Q, n, k) 



e(Q,n) + (HQ - Mvq) — 
m 



Ml M \ n 2 
2 \ ttlU ] m 2 



Substituting this expression into Eq. (|ST)|) we find 

/ Mv 2 

T Q = (d n e Q )d n {v 2 -vl)-{v 2 -vl)[dle Q + — T 

\ n o 



1 

mj 
2d 
no 



y n £Q 



Mv 2 
n 



ft 2 Q 2 d 2 



(TiQ - Muq) (vQd n v - vd n VQ) 



(66) 



This expression can be viewed as a generalization of our 
earlier result (|54[) to the case of impurity with non-zero 
mass M. The full expression for the rate of scattering of 
such an impurity is given by the combination of Eqs. (|66[) 
and ([26]) . The same result can, of course, be obtained by 
generalizing the approach based on Lagrangian variables, 
Sec. 02 to the case of a massive particle. The respective 
calculation is outlined in the Appendix. 

As a simple check we consider an impurity completely 
decoupled from the liquid. In this case the scattering 
probability must vanish, and we expect to obtain Tq = 0. 
This is easily verified by substituting <9„£q = 0, vq = 
HQ/M and m* Q = M into Eq. ([66]). 

The problem of dynamics of a mobile impurity in 
a Galilean invariant fluid was recently addressed by 
Schecter, Gangardt, and Kamenev.— Comparing our ex- 
pression (1551) with the Eq. (75) of Ref. [l6| we find an 

agreement, provided that their matrix element is 

related to Tq as 



n 



T 



Q 



mv {v 2 



,2 ^2 ■ 



(67) 



The authors of Rcf. characterized the impurity by the 
number of particles N it expels from the fluid and the 
supcrfluid phase $. The relation (|67|) was obtained by 
expressing N and <!> in terms of the energy and velocity 
of the impurity with the help of Eq. (21) of Ref. [l6l 



Dissipative dynamics of holes and equilibration 
of one-dimensional quantum liquids 



D. 



In Sees . llTI and Mil we have evaluated the rate Wq^q+sq 
of scattering of a hole excitation by the low energy 
bosons. Our results (f26|) and (|5"Tj) enable one to study the 
dynamics of the hole excitation, provided the occupation 
numbers N q of bosonic states are known. In the most in- 
teresting case when the Luttinger liquid is in equilibrium 



at temperature T, the rate Wq^q+sq falls off exponen- 
tially at 1 8Q | >■ TJ (v — | Vq I ) . Ifa single hole with energy 
e 3> Tv/(v — \vq\) is excited in the Luttinger liquid, the 
collisions with bosons give rise to a gradual change of its 
momentum at the rate 



F 



dt 



HQ 



hSQ W QiQ+S q dSQ. (68) 



Using our most general result (fBTj) for the scattering rate 
we obtain the force acting on the hole in the form 



2vK 2 T 2 Q T i {v 2 



v q) v Q 



15ft 5 



(v 2 — Vq) 5 v 2 



(69) 



The negative sign in this expression indicates that as a 
result of scattering by the bosons the wavenumber of the 
hole approaches one of the Fermi points, Q — or Q = 
2irno, i.e., the hole is eventually absorbed into the boson 
gas. 



As we stated in Sec. IIV CI our results can also be ap- 
plied to a mobile impurity in a Luttinger liquid. An ex- 
pression for the force acting on such an impurity in the 
case of Galilean invariant system was found in Ref. [l6|. 
Using Eqs. (|5"5")) and (16"T|) we find that our result (f6"9")l 
recovers Eq. (73) of Ref. [lg. 

The collisions of the hole with the bosons result in a 
stochastic motion which should be described in terms of 
the hole distribution function Jq. Assuming again that 
the bosons are in thermal equilibrium and Q is sufficiently 
far from the Fermi points, the evolution of the distribu- 
tion function is controlled by the collision integral 



dt 



dQV Q W Q}Q ,-f Q ,W Q ,, Q ] 



(70) 



One can now use our results for Wq.q+$q to study the 
evolution of the distribution function towards the equi- 
librium 4 0) = e~^Q)/ T . 

The collision integral ([70)1 takes a particularly simple 
form for Q in the vicinity of Qq = 7rno, where the en- 
ergy e(Q) takes the maximum value, Fig. [2] In this case 
the typical change of energy in a collision Se ~ vqSQ ~ 
(vq/v)T -C T. Thus the distribution function changes 
very little after each collision. This enables one to bring 
the collision integral (1701) to the Fokker-Planck form 



dt " 



f = -d Q [A(Q)f~-d Q [B(Q)f] 



(71) 



where 



A{Q) = J 6QW Q , Q+SQ dSQ, (72) 

B(Q) = J{8QfW Q , Q+5Q d5Q, (73) 

see, e.g., Ref. H3- 

The Fokker-Planck equation for the evolution of the 
hole distribution function has been applied earlier to the 
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problem of equilibration of a one-dimensional system of 
interacting electrons ^l ~ 19 i 31 At low temperatures the col- 
lisions lead to relatively fast thermalization of excitations 
near each Fermi point. However, the full equilibration of 
the system includes exchange of particles between the 
right- and left-moving branches, which equilibrates the 
respective chemical potentials. This process involves dif- 
fusion of a hole excitation in momentum space from one 
Fermi point to the other In the case of arbitrary in- 
teraction strength, when the one-dimensional system is 
treated as a Luttinger liquid, the equilibration rate is 
expressed asil 



3hnlB (hv 
y / 2nm*T \ T 



-A/T 



(74) 



Here A, m* , and B are given, respectively, by the hole 
energy e(Qo); effective mass rn,Q o , and B(Q ) at the max- 
imum Qq = nn . Our results for Wq^q+sq allow us to 
use Eq. (|73[) to evaluate B for any Luttinger liquid, re- 
gardless of the presence of Galilean invariance. Using 
Eq. ((SJ) we find 



4tt K 2 T 2 
B = — 



Qo 



15 hJv 1 



(75) 



In the case of a Galilean invariant system, Eq. (fT5"j) re- 
covers the result of Ref. Il7l. 



E. Integrable models 

The excitation spectrum of a one-dimensional quan- 
tum liquid can be obtained exactly for integrable mod- 
els. A simple example of such a model is that of spin- 
less fermions of mass m interacting with potential de- 
caying as inverse square of the distance between parti- 
cles, V(r) = g/r 2 . The hole spectrum of this Calogero- 

.32 



Sutherland model is given bj» 

h 2 X 



e(0.») = -z-Q(2nn-Q), 
Zm 



(76) 



where A(A — 1) = gm/h 2 . The velocity of holes is easily 
obtained from Eq. ([76]). 



(77) 



v(Q, n) = —(mi - Q), 
m 



and the velocity of bosonic excitations v = v(0, n) = 
■nh\n/m. Upon substitution of these values into the ex- 
pression (|54p one finds Tq = 0, which points to the ab- 
sence of scattering of holes by bosonic excitations. 

The absence of scattering of excitations is widely 
believed to be a universal property of integrable 
modelsj 12 i 15 ~ 19 i 33 ~ — Indeed, integrability implies the 
presence of a large number of integrals of motion, which 
precludes scattering of excitations. Given the expressions 
(12(31) and (I5T1) for the scattering probability, we expect 
that Tq = for any integrable model. 



V. SUMMARY 

To summarize, in this paper we developed a phe- 
nomenological theory of scattering of hole excitations in 
one-dimensional spinless quantum liquids. We expressed 
the scattering rate in terms of the spectrum of the hole 
and its dependence on the fluid density and velocity. We 
considered liquids which may or may not possess Galilean 
invariance. 

Our approach is based on the concept of Luttinger liq- 
uid and thus applies at any strength of the interactions 
between the particles. We used two alternative descrip- 
tions of the liquid, based on either Lagrangian or Eulerian 
variables. The former has the advantage of simplicity 
when applied to Galilean invariant systems. The latter 
is somewhat more complicated but provides a more nat- 
ural description of systems without Galilean invariance. 

Our most general result for the scattering probability 
of a hole Wq 7 q + sq is given by Eqs. ([ST]) and ([ST)]) . It 
simplifies considerably for Galilean invariant systems. In 
this case Wq^q+sq is given by Eq. ([26]) where Tq may be 
expressed in two equivalent ways. It takes the form (|2"5)) 
if the energy of the hole is given as a function of the La- 
grangian wavenumber. This description arises naturally 
in the limit of strong repulsion between particles, when 
the system forms a Wigner crystaL 18 i 19 Alternatively, if 
the energy of the hole is known as a function of the Eu- 
lerian wavenumber the expression (|54p for Tq is more 
convenient. 

In our theory the hole excitation is treated as a mo- 
bile impurity. Consequently, our result for the scattering 
probability (IBTl) and (|50j) applies not only to the hole, 
but any impurity, including a foreign particle introduced 
into the system. In the Galilean invariant case, dynam- 
ics of such a particle was studied recently by Schecter, 
Gangardt, and KamenevJ^ We have verified that their 
expression for the force acting on the particle agrees with 
our Eq. ([55]). 

The scattering of intrinsic hole excitations controls 
equilibration of one-dimensional quantum liquids. Apart 
from fundamental importance, the equilibration de- 
termines the conductance of long uniform quantum 
wiresj 31 ' 36 Previous calculations of the equilibration 
rat©i£~— assumed Galilean invariance. Our discussion 
in Sec. II V Dl extends these results to the general case. 
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Appendix: Description of a mobile impurity in 
Lagrangian variables 



In Sec. HVCI we discussed scattering of a massive mo- 
bile impurity in a Luttinger liquid. Our approach was 
based on the conventional description of the system in 
terms of Eulerian variables. In this Appendix we de- 
rive the main result (I66[) using Lagrangian variables, see 
Sec. HU We start by generalizing the expression (TTTj) for 
the Hamiltonian of the excitation to account for non- 
vanishing mass M of the mobile impurity. 

Following the approach of Sec. UH we consider an ele- 
ment of the liquid of length Ay which contains the im- 
purity. We write the energy of the element as a sum of 
the kinetic energy of center of mass of the system and 
the energy in the center of mass frame, 



AP = 



(AP) 2 



2(mn Ay + M) 



+ U(n)n Ay+e(P Y /h,n), (A.l) 



cf. Eqs. ((TJ and (JTT])- Expanding Eq. (jA.lj) to first order 
in M I (mnoAy) -C 1, we obtain 



(AP) 2 

AH = -i >— + U(n)n Ay 

ZmnoAy 

+e(P Y /h,n)- \m(-^- 



(A.2) 



Here Py and AP are the momenta conjugated to the La- 
grangian coordinate Y of the impurity and the center of 
mass coordinate R of the element, respectively. Because 
of the non- vanishing mass of the impurity, the latter is 
no longer equivalent to the displacement u of the fluid, 



R 



mn Ay[y + u(y)] + M[Y + u(Y)} 
mnoAy + M 



(A.3) 



Here y is the Lagrangian coordinate of the fluid element 
containing the impurity, y — Ay/2 < Y < y + Ay/2. A 
small change in the positions of the element du and the 
particle dY results in the shift of the center of mass by 



dR = du 



M[l + u'{Y)\ 



dY. 



(A.4) 



mnoAy + M 

Unlike the case of M — considered in SecHJ the center 
of mass position is affected by the motion of the impurity. 
Despite that, the momentum AP is still expressed in 
terms of the operator d u , 



AP = -^<9fi| y=const = ~ihd u = Ayp(y), 



(A.5) 



where p(y) is the momentum of the liquid per unit length 
introduced in Sec. Ill Al On the other hand, the mo- 
mentum Py — ~ihdy\ji =cons t of the relative motion of 
the impurity and the liquid differs from the momentum 
HQ = —ifrdy\ u=const of the impurity in a stationary liq- 
uid, 



Py 



M[l + u'(Y)] 
mnoAy + M 



(-ihd u ) 



(A.6) 



We now neglect in the denominator the mass of the im- 
purity M compared to the much larger mass mnoAy of 
the fluid element and use Eq. (|A.5j) to exclude d u . This 
yields 



M P (Y) 

Py - hQ — , 

m n(Y ) 



(A.7) 



see Eq. Q . Since p(Y) is the momentum density of the 
liquid, at u' = one can interpret the last term in the 
right-hand side as — MV, where V is the velocity of the 
liquid. 

Our theory is constructed in terms of the displacement 
of the fluid u rather than the center of mass coordinate 
R. Thus we substitute Eqs. (|A.5I) and (|A.7[) into the ex- 
pression (|A.2j) for the Hamiltonian of the fluid element. 
Then the first line of Eq. (fP|) gives the energy of the liq- 
uid without the impurity, cf . Eq. , whereas the second 
line gives the Hamiltonian of the impurity in the form 

= ( Q _Mp { Y) ^ * 



H, 



m Hn(Y) 









mno 



(A.. 



At M = this Hamiltonian recovers that of the massless 
hole, Eq. (fTT|). 

To obtain the probability Wq^q+sq of scattering of the 
impurity from state Q to Q + SQ per unit time, we re- 
peat the steps outlined in Sec. Ill Cl using the Hamiltonian 
(|A.8[) instead of (fTTj) . The first step is to perform the uni- 
tary transformation U^Hl+H^U of the form (fT5|) . The 
latter transforms the relevant operators as follows 

tfu'(y)U = u '(y) + f p 6(y-Y), (A.9) 
U^p(y)U = p(y) + f u S(y~Y), (A.10) 

&QU = Q + \fuu'{Y) + ^f pP {Y). (A.ll) 

We then choose the coefficients f u and f p such that the 
linear in the bosonic fields u and p contribution 

pW = (mnav 2 f p + v Q f u - n d n tQ) u(Y) 
f fu , £ Mvq' 



\mno 



p(Y) (A.12) 



to the transformed Hamiltonian U^(Hl +Hi)U vanishes. 
This gives 



VQ(n d n eQ - Mv 2 



fp 



n d n e Q - Mv 2 Q 
mno(v 2 — v 2 ^ 



Q> 
(A.13) 

The next step is to expand the transformed Hamilto- 
nian to second order in bosonic fields, 



3a/ p + n d n e Q + inga 2 e Q 



-VQfu - n d n VQfu 



f 2 

J 11 



2 m 



u' 2 (Y) 



(mnof P -M) 2 M 



2m 



Q 



Q. 

P 2 (Y) 
m 2 n.Q 



.(A.14) 
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Here ellipses denote the omitted terms of the form 
u'(Y)p(Y) which generate coupling of the impurity to 
two bosons on the same branch and are therefore not 
relevant for our scattering problem, see Fig. [5J 

We now substitute the expressions © and © into 
(|A.14[) and extract the matrix element corresponding to 
annihilation of boson q± and creation of boson 172 , 



>K/ 1 9i 92 1 f Q , n 1 2 02 

r — 3a/ p - n d n e Q - -n d n e Q 

rrinoLv \ 2 



+vqU + n d n v Q f u + 



f 2 
2m* r 



(mnafp - M) V Mv 



2m* 



(A.15) 



Substituting relations (p?T) and (|A.13p we bring this ex- 



pression to the form (1221) with 



(d n e Q ) d n (v 2 - 4) - (, 2 - 1$) (flfo + ^) 



1 Mv< 
5T cW 

H VVn(vd n VQ -VnO n V). 

n 



(A.16) 



This expression generalizes our earlier result (|2"5)l to the 
case of a massive mobile impurity. One can now substi- 
tute Eqs. (ITU1) and (|T5|) to express Tq in terms of the 
physical energy e of the impurity in the liquid. This pro- 
cedure transforms the expression (|A.16|) to the form (prJj) . 
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